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ABSTRACT 


Flow visualization of low speed flow has been obtained 
in this thesis by the injection of a thin stream of gas. 

This injected gas has a refractive index which differs from 
that of the uniform flow in the wind tunnel. A frozen fringe 
hologram of .the introduction of this gas indicates regions 

of fringe shift. A pathline is produced by connecting points 
of maximum fringe shift. Pathlines generated in such a man- 
ner around a cylinder agree with the stream functions of 
potential flow theory. 

A rapid means has been developed of data reduction of 
fringe information from a photograph of a reconstructed holo- 
graphic interferogram. Data points are gathered using a 
digitizer and a mini-computer. The array of data points so 
obtained can be examined and corrected if necessary and fin- 
ally expanded by means of Lagrangian interpolation. The 
expanded array is transferred to a fortran program which ac- 
complisbes the integral inversion and subsequently prints 
out the density field. By using this method of data reduc- 


tion and inversion, a large time savings is realized. 
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I. INTRODUCTION 


Thisethesis Consists of two parts. The first part is an 
investigation of low speed flow visualization by means of 
поТосларвте  ппеегГегошесгу The Mack of Significant density 
changes over a profile at low speeds has made interferometric 
investigations D k t. Tlow.difficult. Іп this study, 
flow visualization at low speeds has been obtained by inject- 
ing a thin stream of gas in the wind tunnel. This gas had a 
refractive index which varied from that of the uniform flow 
in the wind tunnel. A frozen fringe hologram of the flow 
field with the introduced gas indicated regions of fringe 
II T Фе равй ле Gh Бе found by connectingypoints о: 
maximum fringe shift.  Pathlines generated in such a manner 
were compared with potential flow theory by superimposing a 
family of stream functions over the experimental pathline. 
This comparison was done with a cylinder profile in a uniform 
steady flow. 

The second part of this thesis concerns itself with a 
rapid means of data reduction of the large amount of fringe 
data that is obtained in holographic interferometry. A dig- 
itizer and a mini-computer were used to effect this rapid 
jl Sr ага reduction the number of§data points 1S ex- 
panded in the mini-computer by uSing Lagrangian interpolation. 
The whole array of points are also automatically transferred 
| Oo C ran дата reduction program. Phys w&ijisting fortran 


program then inverts the integral cquations to determine the 








density field. In the past, fringe data had been collected 
from a hologram by the use of an optical comparator to deter- 
mine the amount of fringe shift along each plane. The data 
was typed onto IBM cards and fed into a fortran program for 
reduction. The present development of the use of a digitizer 
and the ability to automatically feed data into a fortran 
program greatly reduces the time to obtain a density field 


from a hologram. 





Il. FLOW VISUALIZATION 


A. HOLOGRAPHY AND FLOW VISUALIZATION 

Dennis Gabor in 1948 (1) was the first to invent a meth- 
od of wavefront reconstruction called holography. With the 
development in the early 1960's of a coherent source of ra- 
diation, the laser, the field of holography rapidly developed. 

Holography is the process of storing wavefronts ina 
photosensitive material called a hologram. At a later time 
these wavefronts may be reconstructed from the hologram to 
form an image which has the same depth and parallax charac- 
teristics as the real object. This is accomplished by re- 
cording the interference pattern between the light diffracted 
by the object and a reference beam of light. An optical 
arrangement is set up to split the coherent light source into 
a reference beam and a scene or object beam. Within certain 
restrictions the paths of these two beams are arbitrary. 
The photosensitive plate is placed in such a position that 
the reference beam and the diffracted light from the object 
fail upon it. As stated in Ref. 1, each point on the holo- 


gram contains the full information of the object as seen 


wont bhat point.: Full information constitutes intensity and 
phase information. The phase information is derived from 
the interference pattern. When the object is reconstructed, 


each of these points superimposes its own complete picture, 
which bas depth and parallax. The recording of phase infor- 


Naciona the associated dk Ка рала) 27 hat 
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distinguishes holography from photography which only records 
intensities. 

A frozen fringe hologram is one which has been doubly 
exposed. One exposure records the reference situation. The 
second exposure records whatever perturbation has been intro- 
duced. In the present low speed flow visualization experi- 
ment the reference condition of exposure is the steady flow 
of uniform gas over the two-dimensional profile. After the 
reference exposure an optical component is moved slightly to 
provide fringe lines. The flow field is perturbed by the 
introduction of a thin stream of gas into the uniform steady 
flow ahead of the profile. A second exposure is taken on 
the same hologram plate of the perturbed flow. In the re- 
construction of the resulting hologram, the two recorded 
flow fields interfere with each other and show a deviation 
in the superimposed finite fringes. Since both exposures 
utilize the same optical components, the fringe lines are 
lines of constant differences in optical path lengths. Any 
шпретгвесттлом in the optiealzaualiity of the apparatus is thus 
eliminated. 

When a gas of a different density (optical path length) 
is carefully introduced into the flow field, it is swept 
through the wind tunnel and follows a pathline. Holography 
из "Garremely Sensitive to variations in optical path length 
and Dye using holographic techniques Яя line of fringe shifts 
can be seen in the region of the injected gas. By carefully 


Ile mount of sastintroduced in relation to the 


10 


wind tunnel velocity, a small well defined line of fringe 
Shifts can be produced. A pathline can be created by con- 


necting points of maximum fringe Shift. 


B. EXPERIMENTAL APPARATUS | 
The arrangement of the experiment is shown in figure 1(а). 
This type of set-up is known as a modified Z-Type Schlieren. 
A-Korad Q-switched laser (ruby wave length 6943 A) is mounted 
on an optical rail. A small five-milliwatt helium-neon con- 
tinuous gas laser is placed adjacent to the rail to provide 
alignment of the optics. The beam is expanded and projected 
onto a parabolic mirror. This mirror collimates the beam to 
facilitate easier reconstruction. The reference beam travels 
through the beam splitter and is reflected off two mirrors 
before illuminating the holographic plate. The object beam 
is reflected off the beam splitter and then projected through 
the wind tunnel. Just before the tunnel the object beam 
passes through a diffusing plate and a grid. The actual ex- 
posure time with the Q-switched laser was on the order of 
twenty nano-seconds. This stopped all movement during the 
exposure, but did not prevent movement between exposures. 
The tunnel was mounted Separately from all other components 
to minimize unwanted movement caused by vibration. Alignment 
of the optics was done by projecting the small laser beam 
through the ruby laser head and out into the optical system. 
The initial alignment of the small laser with the Q-switched 
Ec 1ccompisshod with burn holes in photographic paper. 


Tit ndue jS explained in Appendix A of Ref. 2. Tringe 
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lines were formed by moving the beam splitter .003 inches be- 
tween exposures. The beam splitter was mounted on a preci- 
sion X-Y translating stage. 

The wind tunnel is 3.5 inches by 3.5 inches and construct- 
ed out of plexiglass (not of optical quality). The speed of 
the tunnel can be varied from zero to twenty-two feet por 
second. A section of honeycomb was placed at the entrance to 
serve as a flow straightener and as a support for the tube 
which injects the gas. 

The tube used for injecting the gas had an outside dia- 
meter of 3/32 inches. It was connected by tubing to a flow 
meter. This meter allowed visual adjustment of the injection 
gas flow rate. 

Two profiles were used, both constructed from plexiglass. 
Pictures of the profiles and their supports are shown in Fig. 
1(e). A cylinder was used having a length of three inches 
and a diameter of one inch. The second profile was a wing 
having à length of three inches, a chord of 1.5 inches and a 
thickness of 3/16 inches. Both profiles were supported by 
two upright pieces of plexiglass. A grid was lightly inked 
onto the supports to provide a reference during reconstruc- 
tion. The ends of the profiles were painted black for better 
КООШО О оТ their outline.— There was a snug fit between 
the pins, which extended from cach end of the profiles, into 
DNCEEPWONcUubports. This fit! allowed the wing to beset at a 
desired angle of attack and remain there. When the profile 


was adjusted between the supports then the assembly was slid 





into the tunnel. The honeycomb was placed at the entrance 
and the injection tube inserted through the honeycomb. 
Reconstruction was done with a fifteen—-milliwatt helium- 


neon laser. The beam was expanded from .04 inches to .25 


inches by the use of two mirrors. A sketch of the recon- 
struction is shown in Fig. 1(d). A shutter was placed be- 
tween the second mirror and the hologram plate holder. In 


order to photograph the real image the camera lens was re- 
moved. A dark cloth was used to cover the laser and prevent 
any extraneous light from overexposing the film. 
Other equipment used in this experiment: 
1. power supply and circulator for Q-switched laser 
2. eye protectors 
3. Agfa 10E75 holographic emuision plate and Polar- 
oid Type 55 film 
4. auto collimator for internal Q-switched laser 


alignment. 


C. EXECUTION AND RECONSTRUCTION 

Before taking a frozen fringe hologram, the alignment of 
the optical system was verified using the small 5-milliwatt 
laser. While this was being done the power supply for the 
Q-switched laser was turned on. The recommended warm-up 
time for the power supply is five minutes. The circulator 
was also turned on to stabilize the laser head and the etalon 
at a low temperature. A small amount of air was blown onto 
Liewevtitcwesuimaces Lo remove dust. The mirror directing 


the alignment beam through tbo ruby was removed and a 2.6mm 


Ша 





transverse mode selector pin hole was placed in its holder. 
The Agfa 10E75 emulsion plate was placed in the plate holder 
which was under a dark cloth. After appropriate warnings to 
other personnel in the area and donning safety glasses, the 
laser was charged. The wind tunnel was set at the desired 
speed and allowed to stabilize. The black cloth covering 
the holder was removed, the laser fired, and the black cloth 
replaced. Then the beam splitter LE moved .003 inches. 
The injection gas was turned on and the flow regulated. The 
capacitor bank for the laser was again charged and, after 
removing the cloth, fired. The laser, wind tunnel and the 
injection gas were turned off while developine the emulsion 
plate. 
Emulsion plate development was standard: 

1. Kodak D-19 developer for five minutes 

о = сег тетаста = фор Bath for thirty seconds 

3. fixer for five minutes. 
The plate was washed in water followed by a rinse in Kodak 
Photoflo wetting agent. A hot air dryer was used to dry the 
plate. Occasionally the plate was overexposed and a bleach 
was used. Whenever this was done an index matching solution 
of xylene was also used to reduce the noise level of the 
hologram. Further information on bleaching and index match- 
mie oS avalltable in section 10 е > ог Бет: 1. 

Reconstruction was done with a fifteen-milliwatt Stabi- 

lite laser using the setup explained in section B and 


ПИЕСИ ИО а) The virtual imago was quickly 
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reconstructed to assess the quality of the hologram. This 
was easily done by using the collimating lens of the Stabi- 
lite laser and looking through the hologram. If the quality 
of the plate warranted a photographic record then the plate 
was reversed in its holder. The collimating lens was re- 
moved from the laser and the beam was slightly expanded to 
reduce speckle. The virtual image was projected onto the 
frosted glass at the back of a lensless camera. For a dis- 
cussion of real and virtual images the reader is directed to 
Ref. 1. An f-16 aperature was used on the shuttering device. 
Black cloth covered the camera, holder, shutter and the la- 
ser. The plate holder was adjusted to present the correct 
picture on the ground glass of the camera. A photographic 
record was made using Polaroid Type 55 positive/negative 


11m. 


D. COMPARISON OF RESULTS WITH POTENTIAL FLOW THEORY 
Two-dimensional potential flow theory was used to give a 
comparison of experimental pathlines obtained by connecting 
points of maximum fringe shift. A discussion of this theory 
is presented in Ref. 4. 
The program listed in Appendix Á calculates stream func- 


ONS Of wnitormyilow past a circular cylinder with circula- 


tion. The stream function y is given as: 
a? y Е r 
UNA TU sc In 
у y 1 21 GÀ 
EIN wt 
where  U -» uniform velocity ; qd DIE Te 
а = cylinder radius; ле C reni nti 


To 





An example of the graphical output of a cylinder with circu- 
lation iS Shown in Fig. 2(c). In the experimentation for 
this thesis there was no circulation around the cylinder. 
When the photograph of the cylinder in Fig. 2(a) is super- 
imposed over the stream functions calculated from the program 
in Appendix A then Fig. 2(d) is the result. Notice the close 
agreement between the stream function labeled ß and the 
points of maximum fringe shift over the profile. Downstream 
of the cylinder the shift in the stream function B is be- 
lieved due to the heavier carbon dioxide used as the injected 
gas. As the wind tunnel speed is increased the point of di- 
vergence of the injected gas with the stream function is 
shifted further downstream. It was found that for very low 
шомло speeds the fall of the carbon dioxide occurred at 
earlier points. At wind tunnel speeds of 22 feet per secon 
the line of fringe shifts caused by injecting carbon dioxide 


agreed favorably with a stream function over a cylinder. 


E. DISCUSSION OF FLOW VISUALIZATION 

Three types of injection gases were investigated. Heliun, 
argon, and carbon dioxide were introduced лп оле Лоу Тїсо1а. 
Helium and carbon dioxide provided an acceptable fringe shift. 
The most advantageous gas was carbon dioxide since its high 
molecular weight prevented the injected stream from rapidly 
(I l jJ i E i Or eje pie ај shows that this*iringe 


shift is small and well defined over the profile. 
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TO accurately define the pathline the injected gas must 
be swept through the tunnel. If the flow rate out of the 
injection tube is too high, then the subsequent path of the 
gas does not follow a pathline. To ensure pathline accuracy 
the velocity of injection gas is matched as closely as pos- 
Sible to the velocity of the tunnel. Both the tunnel veloc- 
ity and the injected gas velocity were 22 feet per second. 
Velocity matching required calibration Of both the wind tun- 
nel and the injection flow metering device. Once the veloc- 
ities are matched then the amount of gas injected is dependent 
on the injection tube diameter. At low speeds a tube which 
had a diameter 5% of the height of the tunnel was able to 
define a fringe shift throughout the hologram. At higher 
tunnel speeds a smaller tube can be used since ine gas fiow 
rate is increased and more gas is introduced to cause a 
shift. Fig. 2(b) is an example where the flow velocities 
were not matched. Notice in the figure how the injected gas 
is blown onto the lower surface of the wing. The Kutta Con- 
dition is easily recognized in the circled area at the 


trailing edge. 


F. CONCLUSION DERIVED FROM FLOW VISUALIZATION 

шие Sensitivity Of nolography to Slight variations#in 
density makes this flow visualization technique ideal for 
low speed flow investigation. Depending on the nature in 
which you inject the gas one can emphasize different features 


им пе low. Although it wet not done for this thesis, 
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perhaps vortices can be visualized by injecting more gas. 
As shown by Fig. 2(a), when the velocity of the injected gas 
is matched to that of the wind tunnel the resultant line of 


fringe shift makes a pathline visible around the profile. 
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III. DATA REDUCTION 


A. RAPID DATA REDUCTION TECHNIQUE 

Due to the large amount of data generated by holographic 
interferometry, a rapid means of data reduction is desired 
in the determination of a density field from holography. To 
date the emphasis has been on developing a better fortran 
program which inverts the integral equations in a calculation 
of a density field. A comparison of the different methods 
used to invert the integral has been discussed by Sweeney 
[Ref. 6]. This paper addresses the problem of data acquisi- 
tion and transfer of these data to a fortran inversion pro- 
gram. A digitizer is used in determining the fringe shift 
in an interferogram obtained by reconstructing the frozen 
fringe hologram. A mini-computer processes the data using 
Lagrangian interpolation, and then stores it. At an appro- 
priate time the data are recalled from the storage tape cas- 
sette in the mini-computer and automatically transferred to 
Ame ree computer for density calculation. fm the previous 
method of data reduction the fringe data had been collected 
by the use of an optical comparator which determined the 
amount oi fringe shift along each plane of a hologram. The 
fringe data were then typed onto IBM cards and submitted to 
Morel tram provram in as laree computer.) ine present develop- 
ment of the use of a digitizer and the ability to store data 
S Automatically franser 1t to а fortran program 


proatly reduces the time needed to obtain a dcusily ficld. 


D 








To obtain a photograph of fringe data, a doubly exposed 
hologram is taken of the aerodynamic phenomenon. The first 
exposure is of the reference condition which usually is a 
no-flow situation. Flow is then introduced into the area in 
question and a second holographic exposure is taken. Finite 
fringe lines are established on the hologram by a slight 
movement between exposures of an optical component. The 
holographic plate is developed and dhën the scene is recon- 
structed using a laser source. A photographic record can be 
taken of the real or the virtual image. For a further ex- 
planation of holographic interferometry the reader is direct- 
ed to Ref. 5. 


The actual integral inversion in the fortran program can 


be done in many ways. A comparison of different methods has 
been discussed by Sweeney [Ref. 6]. Some of these methods 
are: 


1. Representation of the integral by orthogonal poly- 


nomials 
2. Fourier Expansions 
3. Fourier Transforms 
4. Fast Fourier Transforms. 


The overwhelming amount of data generated by holographic 
interferometry requires a rapid method of data acquisition, 


evaluation, storage, and transfer. 
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B. EQUIPMENT DESCRIPTION AND OPERATION 

This data reduction technique uses a Hewlett Packard dig- 
itizer (HP9864A), calculator (HP9830A) and printer (HP9866A). 
A photograph of the apparatus is shown in Fig. 3. An acous- 
tiec coupler was also used to allow the calculator to act as 
an on-line terminal with the Naval Postgraduate School IBM 
360 computer. The power supply for the digitizer was mounted 
under the table and the printer, as can be seen in Fig. 3, 
me restingyonstop of the calculator. 

Data were introduced into this equipment with the use of 
the cursor attached to the digitizer. A button is on the 
cursor which defines the origin on the magnetic table. The 
coordinates of single points are entered by depressing but- 
ton "S". Continuous points can be entered by depressing "С". 
This last feature was not used in the program. 

The programmable calculator uses Basic language and has 
an initial memory capability of 1760 words. Additional тет- 
Sry can be added externally or internally to the calculator. 
A tape cassette with 32,000 word storage capacity can be 
ССРО che calculator. Program or date can be stored on 
or retrieved from this cassette. The entire program is on 
the first portion of this cassette. The remainder of the 
Шарпе са Јаре LS used tor store information for up to nine 
pictures, each picture containing twenty fringe lines and 
twenty planes. By employing a Terminal Read-Only-Memory 
(Term ROM) the calculator becomes an on-line terminal for an 


ШЕ СО computer, While gn the Term ROM mode the calculator 
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retains its basic operating ability. It is possible, al- 
though not advisable, for an operator to execute a terminal 
command and while the IBM operation is being performed, use 
the calculator locally. Operation of the calculator as a 


terminal was in accordance with the normal terminal procedure. 


C. PROCEDURE 

The procedure of using this method of data reduction is 
extremely straightforward. A listing of the programs to 
transfer data is given in Appendix B. Each of the six pro- 
grams is stored in a file on a tape cassette. Each program 
is referred to by the cassette file in which the program is 
stored. The arrangement of the programs were designed to 
allow an operator to collect and expand fringe data for up 
to nine pictures. At a later time, when the terminal mode 
for the IBM 360 is available, then the stored data are sent 
to the large computer for reduction. A flow diagram of 
these operations is shown in Fig. 4. A description of the 
procedure used in the data reduction follows. 

A photograph of a frozen fringe hologram was taped to the 
maple of the digitizer. If the hologram did not have a well 
defined grid than a thin transparent grid was taped over the 
EN osraphocbo provide a reference. The calculator, printer 
and digitizer were turned on. The data reduction tape cas- 
SOLLO was placed in the calculator and the command "LOAD" 
was executed. The program was started at point A on Fig. 


4 by executing the command "RUN". Statement 20 requested 
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the operator to enter the time, day, month and year. These 
numbers were introduced into the calculator by typing them 
onto the display. All numbers were separated with commas. 
The display was executed and the program (HP file 0) contin- 
ued to run. HP file O printed the time, day, month, year, 
and the program organization. At statement 190, HP file 1 
was automatically retrieved from the cassette and loaded in- 
Im he calculator. HP file 1 asked for the orientation of 
the grid. The origin and an arbitrary second point on the 
horizontal axis were entered by depressing the cursor "O" 
(for the origin) and then the cursor "S'", By lifting the 
cursor 1 'inch off of the table the established origin was 
erased and the orientation procedure had to be repeated. 
Whenever the calculator is ready to accept a data point from 
Wine digitizer a soft beep is heard. The calculator asks for 
the number of fringe lines to be used from the picture. Next 
the number of planes in the picture is entered. The program 
w able to accept up to 20 fringe lines and 20 planes. The 
equation of each fringe line is entered by entering two arbi- 
meaty Пара points on the line in the linear region. After 
all lines were entered the program prints the mean and stand- 
ard deviation of the slopes of all lines. Also the differ- 
ences from the mean of each slope are printed. The operator 
makes an evaluation as to the validity of his fringe lines 

by comparing the slopes. If any slope varied greatly from 
the mean then all E e Шаа атте тсевкетса, AD example 


Ne printout 1sfshown in Appendix С. If the fringe lines 
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are satisfactory then 'CONT" was executed and HP file 2 is 
loaded into the calculator. The picture number and the num- 
ber of data points in each ne are entered. As each data 
point is entered the elevation and fringe shift are printed. 
AL Conclusion of each plane the array of elevations and 
shifts are stored on a cassette tape file. The first number 
in the array is the number of data points in that plane. 

The second number in the array is the elevation of the first 
data point. The third number is the fringe shift for the 
first data point. Elevations and shift alternate through 
the remainqer of the array. 

The reference file is filled after acquiring all the data 
points in the picture. The reference file must be compatible 
math the fortran data reduction program. A program based on 
Fourier Transform methods, [Ref. 3], was modified to accept 
data input parameters in a Specific order. After the refer- 
ence file is filled, it is stored in the cassette under the 
appropriate picture number. The parameters in the reference 
file are shown in Appendix B, file 2 starting with statement 
430. The operator was at point B in Fig. 4. Either ће 
can gather more data from different pictures or expand the 
data for fortran reduction. 

HP file 4 is loaded for data expansion. Expansion is 
done using Lagrangian interpolation. All of the data points 
from the plane to be reduced are printed out. The program 
asks for the data point numbers which will be the basis of 


Шар аттау о: expanded data. The number of points to use for 
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curve fitting should be less than seven. If a greater number 
of points are used then the computational time greatly in- 
creases but more importantly the curve fit to the data points 
is erratic. If more than one picture needs to be reduced 
then the same point numbers are automatically used for data 
expansion of each picture. The first number in the expander 
array is the number of data points to follow. The second 
number is the number of the plane in the picture. Before the 
expansion of the array all positions were initialized to 
9.99999. This is shown in File FTO3FOO1 of Appendix E where 
the actual data that has been fed in are indicated and those 
positions in the array that are not expanded data are 9.99999. 
The reason for initializing the array is for convenience of 
data transmission. The absence of the nines following the 
array should be investigated by reading out the stored data. 
After the expanded data points are stored in the cassette 

the operator loads HP file 3. 

File 3 is the program which automatically feeds the data 
into the fortran system. Directions for its use are printed 
by executing "LIST 10,110". To reduce the transmission time 
the expanded data is sent to the IBM 360 in arrays of six 
numbers. The arrays of six numbers are placed on a storage 
. disk in the IBM 360 system under a file labeled "FILE FTO3FOOL". 
шие Тогъгап reduction program calls FILE FTOSFOO1 for its 


data input. 





D. DISCUSSION 

The previously described technique permits extremely 
карта data aCGuiSivion and transier. The ability to rapidly 
acquire, transfer, and store fringe data in large quantities 
allows the operator to process many holograms. To determine 
a density field around an asymmetric object in a wind tunnel, 
as many holographic views as possible are used. Each of 
these holograms is taken at a ee angle and provides 
fringe data as seen from that angle. This program has the 
ability to process nine holograms of 400 data points each. 

An improvement of this data reduction method can be made 
by incorporating a plotter peripheral and by developing a 
better interpolating routine. At present the expanded data 
is based on six fringe shift data points. When more than 
Six points are used the expanded array becomes erratic. One 
would like to utilize as ius of the available 20 data points 
as possible on each plane. A new interpolating routine could 
be developed to accept more than six points and still yield 
accurate results. The expanded array and the original data 
points could be compared on the plotter. From the graphical 
output of the plotter an improved subset of data points could 
be chosen to expand the array. By incorporating this scheme 
not only the best order polynomial is used to expand the data 
but also the best set of data points. 

With the ability to evaluate and store up to nine sets of 
fringe data on a cassette, valuable operator/computer time 18 


ШО Eras terrineg and coriceting fortran inpul data. 
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E. CONCLUSION 

A procedure for rapid fringe data collection, transfer 
and storage has been accomplished. With slight modification, 
this procedure is applicable to all fortran data reduction 
techniques. With one photograph of 20 fringe lines and 20 
points, the use of a digitizer was 20 times faster than that 
of an optical comparator. This was assuming the optical com- 
parator was used in conjunction with a programmed desk calcu- 
lator to calculate fringe lines and fringe shifts. The amount 
of time in the transfer of data to the fortran program was 
15 times less by using this technique. When many pictures 
are taken to completely determine the density field around 
an object, the magnitude of time saved is very large. 

The ability to rapidly store and transfer the large 
amounts of data obtained from a finite frozen fringe holo- 


gram results in a considerable time savings. 
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Figure а Profile and supports inserted in tunnel. 


Notice honeycomb and injection tube. 
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Laser 
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Holosram 
Plate Holder 
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Lensless Camera 


Figure 1(e). Schematic of reconstruction 


of a hologram. 








Figure 2(a). 


Photograph of a reconstructed hologram 
Gis a Cylinder in a uniform flow of 

22 fps. Notice the line of fringe 
shifts over the profile caused by the 


injection of a thin stream of carbon 


dioxide. 
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Photograph of a reconstructed hologram 
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Figure 2(c). 


Graphical output of potential flow program 
given in Appendix A. The stream functions 


are for a uniform velocity of 10 ft/sec 


пт акш ште о ШОШО with a clock- 


Se Girculativon ed IE /sec. 


34 





 Z2-DIHENSIDNRL POTENTIAL РЕЦИ | 
STREAM FUNCTIONS 





Figure 2(d). Fringe lines of Fig. 2(a) superimposed on 
graphical output of stream functions around 


a cylinder. 








Figure 3. 


A view of the Hewlett Packard digitizer, 


printer and calculator. On the digitizing 
table is a finite frozen fringe hologram 
with a transparent grid and the electronic 
cursor. To the rear of the table is the 


acoustic coupler. 
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Figure 4(a). Flow diagram of data transfer program. Dashed 
lines are manual operations. Program file num- 
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Figure 4(b). Flow diagram of the data transfer program. 
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APPENDIX A 
Listing of Potential Flow Program 
This program calculates stream functions around a cylinder 


with circulation. 
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APPENDIX B 


Listing of Calculator Data Transfer Programs 
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tion of the files in the tape cassette. By executing "Run 


200" the stored data of any picture may be printed. 
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For each data point entered this 


The array of fringe shifts from each plane are 
A reference file is also filled and stored 


program calculates the fringe shift from the appropriate 


БР сіл of HP File 2. 
stored on tape. 


fringe line. 
оп tape. 


n =. 
as DIM 
I = 
2 — 
= 1 
— [o 
A .., 
Bansa NOI 
ee 
Du ES 
T eg 
...... “ PI; 
iid === 
UE i ae 
—=. Se 
=== "mimo эш 
s. „==. 
AD с=т 
. - г” - i 
— s ~ " = 
wm š š سے‎ 
c > em 
pU ОО < 
а PES 
"t A 
.... = -—7 
irr an == 
. - ma. 
a : 
m 


wo EI = 
EE == —! 
== O 
nm E рос 111 5: 2 
= = г. СЕЗЕ 
: وا ا‎ i 
: => 
СЕН = Im: = E 
= Lr - F oh 
par Sar = = 1:1 m рог: 
A y =s L: zr o 
و س‎ р So Can en 
ed و‎ : Dar nks: буз 11 т. 
3 — = s 5 ==. Pac 
p — — - cR SUPPE АНИ и SINE 
РБ. Dan ==. E 81 са EE = 
aa r F L=: — Pi ee 
Zea sm Ғы Dorm UM NN I do ud FRE 
mS S: B Doc e eur е 
2227 42722 ЧЕ tra Ееее а Е T gem 
л Dr) hen: ee by ee i e 
کت‎ 2 EET ou жиы шша ui am =. 
pu Кез == = P Ж Ж = — 
p HR = ij. pet - г cov E 
— — م‎ — — > — — == а= م‎ = 
аа s - Li E a EM 
и — fe : mm i 
z š an ale nen = ir == г-- — 113: oa 
ms n "im — Fe uum 
s š ізі pU PS . M E — Е 
m MEE t BE Бада Ho fe d. SITE арт 
afk ESS н Lil ыы JE Ё TU rf. с E 
D u - т-- Е rn =, = TT p. а: 
LL on ша са ESI n C n Bu СЕРЕ: 
rn == —. ru taf = = 
SINE k oo ca Lu L; wis P "n 
PORE nt EE 1 mn лт" НЕГЕ: - = 
i-i mala na Y E. EE fe IE L.L. EN = rS 
S mem oi 5] — E ГІН a: en жи mu Lim 
ғ = اس‎ "E i nd E =" Біз ‚кы nn т» Eb Е 
тас 2-2 - 22 cd us = + He þe- E دہ‎ 
س م د ا‎ T Ты” i НЕ fF ет Ze 
IH И 3 in DN +: „m T qv} a ri 
beim + `"; —, pa PL OUS AX —. — — = e mm m 
mb uo S с - ===: E во во i 4 == rit t РЕ 
on) as Б ты" سس‎ ТЕ E пра | „L: si- H E qur = Bue ae! 
سا‎ ізіме 2-44 Lr. T 22575 i MINES °° = i Te og“ HEHE ИЕ 
THE "E iri E г: э ШЕРІ ЕЕ: =. le 7 1x3 n 1 
} sl pmi $i r 1: Beit ж. = _ т Ta 
— E и > "2 2 .-- е +. Б = = хе", 
pori; IH Ma ЕА = + Li, I, LI r = r Ts VT: : — =, — P 
u — к=з = = д mpm. m9 w. ES Е сыр a س‎ 
E 2 та: DOE TUS s s s. и aT . 
іліге: TUN EN На 1 í да Ee RE i = ры at Li EE E EE 
Š 2 y — iiia - صم س‎ ë انت‎ — 5 — — я 
= = oi int LEIJ - вен oo TEE GG = = = کن‎ i - HS - = Е Saat 
- Ри = з т - = - gera Prior с Mo li з і а ает р = Less 11 به‎ 
= Lim i с x m ER Ls Ld ci т ii ai 1 рея 3 13 т rf 
š я i : a Е D = Zo = == Б = 
سم ست‎ — == I A = - 2; ES dp + Ho 4 Мә ote 1 
: ї i Б 1 Г 1 НА y Li! en er Мр | ВЕ un ٣ 
—— — = E г. Е : سے و‎ == Е Е : n p 
Е Е ул E ME 1 = тер ГҮҮ ==: Ez 1 " ът = = = LL 
аа -- кс к Pin i асыш шс n Е E: H d 4 Ва А 4 аты mc erue uec 
5 SS к = ur E от. : гал: у + M c e : 1 11 ; š+ A 77: 
мк тан аре ка а m u B. vor OE MA Li gH H4 rs 
u. anl E DECUS ` Вере, == : و س ي‎ t =... = === 55 == | . > = EN а t = 
YET ши Е EIN. c IET ПД a oop SU o8 g dE i: E >. e: i ea Е = == в 
NIME I ELE шы ss ús; r s Li. LD ID D S C и sI S د ا‎ JL 5 E VC Gp ee A и Бли 
= = эз ы e = cB 2 = - = — % ті == 
1 DELIA S i TET E = А m į 8 RB 
Di Мы Г — ге E г. LL. — CX LL. LL. B.M рез aou 
ЕЛЕН Ени ен Ена = с HA EA RA SE Sh E 
ii OS EET TE CE EEE P r 5 = 
Т". єз. QUU. PSU. SET. Бул. Qm. s>, EM, oe. iQ. poc $T. =. ,سے کد لل س‎ E == L ipm, p .” - =. г. mm, mT, У 
Toe шы рер шат ІТ ae ee ee OA OSES OTT. 
85 Be He el armen -= MEC E SS = г. "тє э”. T SUM Б и s : - š - 
il a A =p ua hii tir ite gu P Pee ¡e E LE Liu is aio lo mi aio ml a زس ب تسو پس اس‎ io a-i a-i aei ei 





5 LA 


HAIT dag 

ОЕ ПНА sa 

IF Z1-2 THEH zi5 

= 1=21+] 

COTO Jn 

ПАША EE А Досе 


I 


HALT tam 
IN SEO ү IHT. EN 

4 STOP 

CHEM ERASE Cia OES HOHE 


= = 
T arent 
aS 

D 


сан GOTO 395 


Sena а = - - 
y F6 Ta m" А 
Ге Г ағалар іа 


zb dna WE Ун pe 


— 


IH HE rt за 1 


— 
— 


IRE Pt 
КП: 


ей г 


ق 


ب 
— 


e 


org eru 
hn Net fas 


[ua fsa ee — + 


? £3 Y,.1 


] Ea! 
Бал 
PA Па КП: 


г; 
س‎ “== 
— 


r. 


Kern 


А p- г. 
=h mmi тісі. 


Гал 


г. 


План тен 


- ; r- 
гал {эз 
{ы крт 


E EDEN EEES E 
Ai FFF Cot IT 1! 


"Па 


1 


E 
— 


> 


рицар зае 


ma 


Cagnat 


1 UH 
ЛАНЫ 


dite 
t 


E rel гл r 1 fL о.м i 
"ESI 


= 
—- 


r q r ri r rt r 
1 


Ree 8 -— = -- 1 - 
г . =, 
-- 


а Ка 


17 ~ 
и И 
– 


F a TT 
ке Ек Life HUMBER” s 
ЕНЕШІ І1 

Be че е, 

ү cnc c and 

ISE FTF 

ЕККЕН 

= +] ] 

A 

RzETÍ 

Е (j snap Il PT l FD | 1 
R CORRELL TION FLAHE ZI" 


J | 
Б Tio s (jJ ce ET nose Е 
ПЕЕ Еи „СОНТ 


INS FI І 
WATT see 

LUTO ы 

ИТТЕР НЕК ЕР БЕГИНЕ РТ етта 
ПІМ SILS1 TE ЩЧ 

E s 


A DES 


111 


DISP | 

ІНРОТ 

П] ор 
IHP و‎ 


| 
HUNGER OF ТЕНЕТ Тя 


з 


1 


ET 


STIME TE T А и 


дылы ра 15 
imam П 


a 
E 


EI: 
Ее 
а 
aL: 
HT 
Е 
a HILF ZEIT ET FRE а 
EU s 


= i 
i 
Tn ag 
£n NS же. 


DHISEN 
lir TI RT ca 

ШЕН i | P i  LEEHE PUP S] И aa 
J | 

I O l p i e [ib i sih Hii His s 

Bi MNT E c 

ШІ ЗЕН ПЕЕ ЕМЕС а 

ШЫРШЫ ВЕЗИ 

БЕЕК ЕЕЕ s 

INPUT Thad 

Dre eee a S l AHGLE IH MER s 
nbn 
I Ri pu! 
l еј 


STA SYMMETRY Ca" 


46 





=.) =.) 


nr r1 
ым КП: 


zr 
тү 

= 

—L 


--- 
— 


> E 
a 


= 
гае 
D 


1 


yp ply Ey eye 
2 


r: Pa t ERA: 


---. өлі 
ni - 
$ т." 


I 


E etsi + 1 
Ead 7 us es 


MUFE wale 223 7 

А О Hii TO EELUCE THiS PICTURE TMEM...LümHb 5 
БИЕТІ A me DATA FUE HEST PICTURE 
БЕТ" E 

IE AS рез пу ен зи TRE: E 

5 ТИР 

ЕНИ 


47 


E 





This program recalls the array of 


automatically sends the data to FILE 


expanded data from the tape cassette and then the calculator, 
FTOSF001 in the IBM computer. 
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B.5 Listing of HP File 4. This program expands the fringe 
shift data by using Lagrangian interpolation. The array of 


expanded fringe data is stored on the tape cassette. 
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B.6 Listing of HP File 5. The beginning of this program 
compares two pictures at one plane. The second section of 


this program redefines all of the data files. 
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APPENDIX C 


Printed Results of the Calculator Programs 


kest Hron File 0, File 1, and File 2. 
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C.2 Printed output of readout of stored data. 
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APPENDIX D 


HorunsndMoubrout:ne whroh Inputs Data 


E u-EcumPcouPinegwrseun the Fortran data reduction program. 


The title, reference numbers, and fringe shift data are read 


a LE PTPOSE0011: 


1000 
1001 
B9 Z 
T005 
980 
SM 
Bez 
IS 
994 


505 


CUBPOUTINE DATACIT,OA,QB, OS) 

КЕ THIS INPUTS DATA FROM FILE FTOQ3F??? **x**** 
REAL IT(12),QB(6),QS(150) 

INTEGER ОА(5) 


HEAD OS 298 SCIT -1,12) 


Di 
WRITEC6,1000)(ITCID2, 121,12) 
РАВНА ОАО, 1-1,2) 
VEHI, TUI CONGIOS TS 1459 
ADS SUBIO ПЕТ, 6? 
Weve (Cs, 1002) (05(1),1=1,6) 


[D=6*(QA(3)/6+1) 

DO 980 1-1,10,5 

ВЕДО (3,984) (05(1+49-1),9:1,6) 
WRITE(6,1003)(0S(1*J7-1), 071,6) 
FORMAT(1X,6A8) 
FORMAT(1X,518) 
FORMAT(1X,6F10.5) 
FORMAT(1X,6F11.5) 

CONTINUE 

FORMAT (6A8) 

FORMAT (518) 

FORMAT(6F10.5) 
FORMAT(6F11.5) 

1X=9A(3) 

0439525 Wal, Ix 

SOO) 

CONTINUE 

RETURN 

END 
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